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Abstract 

Chemical reaction networks taken with mass-action kinetics are dynamical systems that 
arise in chemical engineering and systems biology. In general, determining whether a chem- 
ical reaction network admits multiple steady states is difficult, as this requires determining 
existence of multiple positive solutions to a large system of polynomials with unknown co- 
efficients. However, in certain cases, various easy criteria can be applied. One such test is 
the Jacobian Criterion, due to Craciun and Feinberg, which gives sufficient conditions for 
ruling out the possibility of multiple steady states. A chemical reaction network is said to 
pass the Jacobian Criterion if all terms in the determinant expansion of its parametrized 
Jacobian matrix have the same sign. In this article, we present a procedure which simplifies 
the application of the Jacobian Criterion, and as a result, we identify a new class of networks 
for which multiple steady states is precluded: those in which all chemical species have total 
molecularity of at most two. The total molecularity of a species refers to the sum of all of 
its stoichiometric coefficients in the network. We illustrate our results by examining enzyme 
catalysis networks. 

Keywords: chemical reaction networks, mass-action kinetics, multiple steady states, 
Jacobian Criterion, injectivity, species-reaction graph 

1 Introduction 

This article concerns an important class of dynamical systems arising in chemical engineering 
and systems biology, namely, chemical reaction networks taken with mass-action kinetics. As 
bistable chemical systems are thought to be the underpinnings of biochemical switches, a key 
question is to determine which systems admit multiple steady states. However, the systems of 
interest are typically high-dimensional, nonlinear, and parametrized by many unknown reaction 
rate constants. Therefore, determining whether a chemical reaction network admits multiple 
steady states is difficult: it requires determining existence of multiple positive solutions to a large 
system of polynomials with unknown coefficients. However, various criteria have been developed 
that often can answer this question more easily. For example, the Deficiency and Advanced 
Deficiency Theories developed by Ellison, Feinberg, Horn, and Jackson in many cases can affirm 
that a network admits multiple steady states or can rule out the possibility [12^ [T^ \T5\ [T9] . 
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These results have been implemented in the CRNT Toolbox, computer software developed 
by Feinberg and improved by Ellison and Ji, which is freely available [13]. Related software 
programs include BioNetX \22\ [25] and Chemical Reaction Network Software for Mathematica 
|21j . For systems for which the software approaches are inconclusive, Conradi et al. advocate a 
two-pronged approach: they first determine whether certain subnetworks admit multiple steady 
states, and, if so, test whether these instances can be lifted to the original network [5]. Another 
criterion for lifting multistationarity from subnetworks appears in [2^. For systems whose 
steady state locus is defined by binomials, a criterion for multistationarity is given by Perez 
Millan et al. in [24] . Finally, a degree theory method which can establish multiple steady states 
can be found in work of Craciun, Helton, and Williams [9]. 

A test for answering the question of multistationarity, which is implemented in some of the 
above-mentioned software programs, is the Jacobian Criterion, which is the focus of this article. 
The criterion, which is due to Craciun and Feinberg, gives sufficient conditions for precluding 
multiple steady states O [Gj [T] [8] . This work has been extended by Banaji and Craciun [H [2] 
and by Feliu and Wiuf [16] . A chemical reaction network is said to pass the Jacobian Criterion 
(or to be 'injective') if all terms in the determinant expansion of its parametrized Jacobian 
matrix have the same sign (see Lemma l3.4p . We prove that for a given network the Jacobian 
Criterion can be checked by examining certain subnetworks (Theorem 14. 9p . 

Our main result (Algorithm 16. ip is a procedure which simplifies the application of the 
criterion in two steps: first, the network is reduced to a simpler network, and then the number 
of subnetworks that must be examined is decreased. As an important consequence, we identify 
a new class of networks for which multiple steady states is precluded: those networks in which 
each chemical species has total molecularity of less than or equal to two (Theorem 15. ip . The 
total molecularity of a species refers to how many non-fiow reactions it appears in, where 
reversible reactions are counted only once and each occurrence of the species is counted with its 
stoichiometric coefficient; see Definition 12.21 More generally. Algorithm 16.11 allows us in certain 
cases to rule out the possibility of multiple steady states simply by inspecting the network. In 
addition, the reduction in the number of subnetworks that must be examined lends support 
to the observation that typically very few terms (if any) in the Jacobian determinant have 
anomalous sign [5]; this topic is explored by Helton, Klep, and Gomez [18] and Craciun, Helton, 
and Williams [9]. A proposed application of our work is to the problem of enumerating networks 
that admit multiple steady states: our algorithm quickly removes from consideration many 
small networks that are computationally expensive to enumerate (see Remark 16. 4p . Finally, 
we provide examples which demonstrate that our approach is complementary to that of the 
Species- Reaction Graph Theorem of Craciun and Feinberg [71 llOj. 

We note that our results are stated for continuous-flow stirred-tank reactors (CFSTRs), 
chemical reaction networks in which all chemical species are removed at rates proportional to 
their concentrations (see Definition 12. 2p . However, even for more general networks in which 
not all species take part in the outflows, our results often nevertheless apply in one of the 
following two ways. First, one can augment a network by any missing outflows so that it 
becomes a CFSTR, and then apply our results to determine whether it passes the Jacobian 
Criterion. Then, under mild assumptions on the original network (such as weak-reversibility; 
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see Definition 12. 2p , the recent work of Craciun and Feinberg allows us to conclude that the 
original network also precludes multiple steady states HI Corollary 8.3]. A second approach is 
due to work of Feliu and Wiuf. Their extended definition of the Jacobian Criterion applies to all 
networks [16], thereby bypassing the need to augment a network by missing outflows. Finally, 
we emphasize that the Jacobian Criterion applies for more general dynamics than mass-action 
kinetics [H [9]; our results also may be used in those settings. 

This article is organized as follows. Section [2] provides an introduction to chemical reaction 
systems. Section [3] describes the Jacobian Criterion, and Section H] lays the groundwork for 
simplifying its implementation. The main result of Section [5] (Theorem 15. ip precludes multi- 
stationarity from any network in which each chemical species has total molecularity at most 
two. Section [6] presents Algorithm 16. H our procedure which simplifies the application of the 
Jacobian Criterion. Finally, Section [7] contains examples such as enzyme catalysis networks 
which illustrate our main results. 

2 Chemical reaction network theory 

In this section we review the standard notation and define 'total molecularity.' An example of 
a chemical reaction is denoted by the following: 

2Xi + X2 ^ Xs . 

The Xi are called chemical species, and 2Xi + X2 and X3 are called chemical complexes. 
Assigning the reactant complex 2Xi + X2 to the vector y = (2, 1,0) and the product complex 
X3 to the vector y' = (0, 0, 1), we rewrite the reaction as y — ?> y' . In general we let s denote the 
total number of species Xi, and we consider a set of m reactions, each denoted by 

for k G {l,2,...,m}, and y/c,?/^ G Z>q, with y^ / y^. We write each such vector y^ as 
Vk = {yki,yk2, ■ ■ ■ ,yks) & ^>0' ^^"^ ^^^^ yki the stoichiometric coefficient of species Xi in 
complex y^. For ease of notation, when there is no need for enumeration we will drop the 
subscript k from the notation for the complexes and reactions. 

A self-catalyzing reaction will refer to a reaction, such as, A + B ^ 2A + C, in which there 
is a species (here. A) that is a self-catalyst: it appears with a higher stoichiometric coefficient 
in the product complex than in the reactant complex. 

Definition 2.1. Let S = {Xi}, C = {y}, and IZ = {y ^ y'} denote finite sets of species, 
complexes, and reactions, respectively. The triple {S,C, TZ} is called a chemical reaction network 
if it satisfies the following: 

1. for each complex y £ C, there exists a reaction in TZ for which y is the reactant complex 
or y is the product complex, and 

2. for each species Xi G S, there exists a complex y G C that contains Xi. 
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We say that a species Xi is a reactant species (respectively, product species) if it appears in 
a reactant (respectively, product) complex of some reaction. A subset of the reactions IZ' <ZTZ 
defines the subnetwork ,C|7^/,7^'}, where C\'jij denotes the set of complexes that appear 

in the reactions TZ' , and S\c^, denotes the set of species that appear in those complexes. A 
network decouples if there exist nonempty sets of reactions IZ' and IZ" that contain disjoint sets 
of species and such that TZ = TZ' (j TZ" . Otherwise, we say that the network is coupled. 

Definition 2.2. 1. A chemical reaction network is weakly-reversible if for each reaction 
y ^ y' G TZ, there exists a sequence of reactions in TZ from y' to y (that is, reactions 
2/' 2/1 2/2 ■■■—>■ yA; yj- ^ reaction y ^ y' & TZ is reversible if its reverse 
reaction y' y is also a reaction in TZ. 

2. A flow reaction contains only one molecule; such a reaction is either an inflow reaction 
— 7- Xj or an outflow reaction Xi ^ 0. A non-flow reaction is any reaction that is not a 
flow reaction. 

3. A chemical reaction network is a continuous-flow stirred-tank reactor network (CFSTR 
network) if it contains all outflow reactions Xi ^ (for all Xi G 5j. 

4. Consider a chemical reaction network whose set of non-flow reactions is given by: 

2/1 -^y'l 2/2^2/2 • • • yi^ y'l 

yi+1 ^ y'l+i yi+2 ^ y'1+2 ■ ■ ■ yi+k ^ y'l+k , 

where none of the first I reactions is reversible. The total molecularity of a species Xi 
in the network, denoted by TM{Xi), is the following integer: 

l+k 

TM{Xi) = J2iy^^+y'n) ' 

i=i 

where yji is the stoichiometric coefficient of species Xi in the complex yj. 

Note that Craciun and Feinberg use the term "feed reactions" for inflow reactions and "true 
reactions" for non-flow reactions. In chemical engineering, a CFSTR refers to a tank in which 
reactions occur. An inflow reaction represents the flow of species (at a constant rate) into the 
tank in which the non-flow reactions take place, and an outflow reaction represent the removal 
of a species (at rate proportional to its concentration). 

2.1 Dynamics 

The concentration vector x{t) = (^xi{t), X2(t), . . . ,2;|5|(t)) will track the concentration Xi{t) of 
the species Xi at time t. For a chemical reaction network {S,C, TZ = {yk — > 2/1-}}) mass-action 
kinetics defines the following system of ordinary differential equations: 

\n\ 

x{t) = 5^Kfcx(t)^H2/fe-2/fc) =: f{x{t)), (1) 
k=l 
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where the last equahty is a definition. Here, we make use of multi-index notation: x(t)^''' := 
xi{t)y''^X2{t)y''^ ■ ■ ■ Xs{t)y'''' . (By convention, := 1). Also, the k^'s denote unknown (positive) 
reaction rate constants. A concentration vector x G M?^q is a (positive) steady state of the 
mass- action system ([T|) if /(x) = 0. 

Remark 2.3. We can scale each concentration Xi so that each outflow reaction (a reaction 
Uk ~^ y'k ^^^^ Uk = o,nd y'j^ = 0) has outflow rate equal to one (Kk = 1 for such a reaction) 
13 Remark 2.9]. We will make this assumption beginning in the next example. 

In the following example and all others in this article, we will label species by distinct letters 
such as ^, -B, . . . rather than Xi , X2, .... 

Example 2.4. Consider the following CFSTR: 

k2 = l fc4 = l 

The total molecularity of species A and B are three and one, respectively. The mass-action 
differential equations ([I]) for this network are the following: 

ki- XA + hxAXB 
h- XB - h^AXB ■ 

Definition 2.5. A CFSTR is said to admit multiple steady states or is multistationary if there 
exist reaction rate constants Kj G M>o whose resulting system ([TJ admits two or more positive 
steady states. 

2.2 Square networks 

Throughout this article, we will be interested in subnetworks in which the number of species 
and the number of reactions are equal. We will call such subnetworks 'square.' 

Definition 2.6. 1. A chemical reaction network is said to be n-square if it contains exactly 
n species and exactly n directed reactions. A network is said to be square if it is n-square 
for some n € Z>o. (A 0-square network is the empty network.) 

2. For an n-square chemical reaction network N = {5,C,7^} with reaction set TZ = {yi — > 
y[ \ i = 1,2, ... ,n}, the reactant matrix and reaction matrix are the integer (n x n)- 
matrices 



dt 
d 

It 



XA 



XB 



yi 
y2 



yi 



and Rjy 



y'l 
y2 - 2/2 



, respectively. 
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3. For a square chemical reaction network N, we define its orientation Or(A^) € { — 1,0,1} 
as follows: 

Or{N) := sign(det(Mjv)det(i2Ar)) . 

As defined, our reactant and reaction matrices are the transposes of tliose of Craciun and 
Feinberg [5]. However, the orientation of a network is unaffected by transposing either matrix. 
Note that a reaction yi — )• y[ is encoded by the vector yi — y[ in the reaction matrix R^, 
which is the negative of its reaction vector y[ — yi which appears in the mass-action differential 
equations ([1]). 

3 The Jacobian Criterion 

The Jacobian Criterion concerns the determinant expansion of the Jacobian matrix of the 
differential equations ([T]) of a chemical reaction network. 

Definition 3.1. Consider a chemical reaction network <8 with s species. The Jacobian matrix 
associated to © is the s x s-matrix whose {i,j)-th entry is the partial derivative -^fi{x), where 
fi is the i-th differential equation in ([T]). The Jacobian determinant o/(J5, denoted by Jac(©), 
refers to the determinant of the negative of the Jacobian matrix, which is a polynomial in the 
variables Xi (for i = 1,2, . . . s) and Ki. 

Note that Craciun and Feinberg use the notation det k)^ for Jac(0)) [5]; here p(g 

denotes the negative of the mass-action differential equations for the network obtained from <d 
by removing all inflow reactions. 

Example 3.2. We now return to the CFSTR examined in Example \2.4\ which we denote by 
&. Its Jacobian matrix is: 

/ -1 + k^XB kr^XA 
V -hxB -1 - k^XA 

Therefore, the Jacobian determinant of © is the following polynomial: 

Ja.(e) = det(l;*;;- = + (2) 

The following result, which is due to Craciun and Feinberg [5], states that the determinant 
expansion of the Jacobian matrix has at most one term in the case of a square chemical reaction 
network. 

Lemma 3.3 (Theorem 3.2 of [5]). Let N be a chemical reaction network that contains s species 
and s directed reactions. Then the Jacobian determinant of N , Jac(A^), is either zero or a 
monomial (in the variables Xi and Ki) whose coefficient is given by the following product: 

det(MAr) • det{RN) , 

where is the reactant matrix and Rn is the reaction matrix. 
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The following result is due to Craciun and Feinberg. 

Lemma 3.4 (Theorems 3.1 and 3.3 [5J). Let & = be a chemical reaction network 

that contains s species, r non-flow reactions, Tq outflow reactions, and ri inflow reactions. Let 
Jac(©) denote the Jacobian determinant associated with (5, which is a polynomial in variables 
xi,X2, ■ ■ ■ ,Xs, Ki, K,2, ■ ■ ■ , Hr- Let Gi, G2, ■ ■ ■ , G ^r+ro+r^-j enumerate the subnetworks of that 

contain exactly s reactions, such that the first C^'") do not contain inflow reactions. 

1. Each term in the expansion o/Jac((5) is the Jacobian determinant of a unique subnetwork 
Gi. Conversely, the Jacobian determinant of each subnetwork Gi is either zero or is a 
monomial in the Jacobian determinant of<3. Ln other words, 

Jac(6) = Yl Jac(G,) = ^ Jac(G,) . 
1=1 1=1 

2. Moreover, if & is a CFSTR, then the following are equivalent: 

• The polynomial Jac(C5) is positive for all x G M^q and Ki > 0. 

• Each term Jac(Gi) is either zero or a positive monomial. 

• Each subnetwork Gi has non-negative orientation: Or(Gj) > 0. 

// these equivalent conditions hold, then (5 does not admit multiple steady states. 

One of the main results in the next section is Theorem 14. 9| which will extend the equiv- 
alent conditions that appear in part 2 of Lemma 13.41 so that instead of checking whether the 
subnetworks Gi have non-negative orientation, we need only check the orientations of certain 
smaller 'embedded' networks (Definition 14. 6p . If any of these conditions hold for a CFSTR 0, 
for instance if Jac(0) has only positive terms, then we say that (S passes the Jacobian Criterion 
(or (5 is 'injective'). 

Remark 3.5. Jac((5) always has at least one positive term when (5 contains all outflow reac- 
tions (such as for a CFSTR). This arises because one of the subnetworks Gi consists of the s 
outflow reactions. Clearly, the reactant and reaction matrices of this subnetwork are both the 
identity matrix, so its orientation is indeed positive. 

The aim of this article is to simplify the computation that must be performed in deter- 
mining whether a CFSTR passes the Jacobian Criterion, that is, whether all other terms in 
the determinant expansion are positive. In fact, there are typically few negative terms [5]; this 
topic is explored in [9l [18] . 

Notation 3.6. Although Feliu and Wiuf have extended the definition of the Jacobian Criterion 
to apply also to non-CFSTR networks 116^ . the main focus of this article is the case when the 
chemical reaction network (3 is a CFSTR. The subnetwork of (5 containing all the non-flow 
reactions but none of the flow reactions will be called the non-flow subnetwork of © and denoted 
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by G. Conversely, any reaction network N = {S^C^TZ} is contained in a CFSTR obtained by 
including all outflow reactions: we denote this CFSTR by 

N:={s, cusyj {0}, nyj{x,^ ^}x^es} ■ 

Example 3.7. We now return to the CFSTR examined in Examples \2.4\ and \3.SX which has 
s = 2 species and r = 1 non-flow reaction. Rs non-flow subnetwork is: 

G := {A + B^2A}. 

The C+Z") = Q = 3 2-square subnetworks of G = (5 = {A ^ B ^ 0, A + B 2A} 
that do not contain inflow reactions are the following: 

Gi = {A ^ 0, S ^ 0} , 

G2 = {A^{), A + B ^ 2A}, and 

G3 = {B A + B ^2A} . 

As noted earlier, the subnetwork comprised of all outflow reactions, Gi, has positive orientation: 
Or(Gi) = 1. Now we examine the reactant and reaction matrices of the remaining subnetworks 
G2 and G3: 

Mg. = ( I J ) = ( _\ J ) (3) 

Mo, - (; ;) no, - 

We see that the orientations of the two networks are Or(G2) = sign (det(MG'2) det(-RG'2)) = 1 
and Or(G3) = sign (det(MG'3) det(i?G3)) = ~1- Therefore, Lemma \3.4\ implies that the CFSTR 
& = G fails the Jacobian Criterion. This is consistent with the computation of the Jacobian 
determinant ([2]) in Example \3.SX which we saw contains one negative term. 

3.1 Reactant and reaction rank 

The Jacobian Criterion requires that each s-square subnetwork Gi of © has either a singular 
reactant matrix, a singular reaction matrix, or both matrices have nonzero determinants with 
the same sign. We now give a name to these cases: 

Definition 3.8. A square chemical reaction network N is said to have full reactant-rank if 
det(Af7v) 7^ 0, and is said to have full reaction-rank if det{RN) ^ 0. 

If a square network has both full reactant-rank and full reaction-rank, then it clearly has 
nonzero orientation. 

Definition 3.9. In a square network N , we say that species Xi is changed in reaction j if 
{RN)j^ ^ 0. 
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In other words, species Xi is changed in reaction j if the concentration of species Xi is 
changed when reaction j takes place. For example, species A is changed in the reaction A — ?> 
B + C and in the reaction A + B ^ 2A + B. On the other hand, species B and C are changed 
only in the first reaction. 

Lemma 3.10. If N is a square network with full reactant-rank, then: 

1. Each species of N is a reactant species. 

2. No reaction of N is an inflow reaction or a generalized inflow reaction of the form — > 
YyiO-iXi (where ai > 0). 

3. No two reactions of N share the same reactant complex. 
If N is an n-square network with full reaction-rank, then: 

1. Each species of N is changed in some reaction. 

2. N does not contain both a reaction and its reverse. 

Proof. This result follows because nonsingular square matrices can not have a zero column, a 
zero row, two identical rows, or a row and its negative. □ 

3.2 Known tests for passing the Jacobian Criterion 

Known sufficient conditions on the reaction network for passing the Jacobian Criterion arise 
from the following results: the SCL Graph Theorem of Schlosser and Feinberg [25], the Species- 
Reaction Graph Theorem of Craciun and Feinberg [71 [8] and the related work of Banaji and 
Craciun [il [2], and results of Helton, Klep, and Gomez |18j . In this article, we present a 
complementary test (Algorithm 16. ip . 

4 Analysis of subnetworks 

Theorems 14.91 and 14.111 in this subsection form the basis of the procedure (Algorithm 16. ip that 
we will introduce to simplify the implementation of the Jacobian Criterion. The following 
example motivates our results. 

Example 4.1. We return to the network G in Example \3. 1I\ Note that in both the reactant and 
product matrices ^ of the subnetwork G2, the first row is the vector (1,0). So, the orientation 
of G2 depends only on the entries in the lower-right of the two matrices. In other words, by 
letting N2 denote the 1-square network obtained from G2 by 'removing' the outflow reaction 
j4 — 7> and 'removing' the species A from the non-flow reaction A -\- B ^ 2 A: 

N2 = {5^0}, 
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it follows that Or(G2) = Or(iV2). Similarly, for the network obtained from G3 by removing 
B and the species B from the non-flow reaction: 

A^s = {A^ 2A} , 

we have Or(G3) = Or(iV3). 

In this section, we make the simphfications performed in Example 14.11 formal : the first rows 
in the matrices Mq^ and (© will be called 'pivotal rows' (Definition 14. 2|) . Theorem 14 . 91 will 
permit us to restrict our attention to the square 'embedded' networks Ni rather than the larger 
s-square subnetworks Gi, and Theorem 14.111 will allow us to remove the species B from the 
network. In general the networks Ni will be smaller and fewer in number, thereby simplifying 
the Jacobian Criterion. 

4.1 Analysis of row and column pivots 

Definition 4.2. 1. For any matrix A, if there exists a row i with one positive entry, Ay > 
(respectively, < 0) for some j, and all other entries zero, Ai^ = for k ^ j, we say that 
Aij is a positive (respectively, negative) row pivot and the i-th row is a pivotal row. 

2. For any matrix A, if Aij > (respectively, < 0) for some i and A^j = for all k ^ i, we 
say that Aij is a positive (respectively, negative) column pivot and the j-th column is a 
pivotal column. 

Notation 4.3. For an s x s-matrix A and indices 1 < i,j < s, we define the associated matrix 
A^^ which is obtained from A by removing the i-th row and the j-th column. 

In the following theorem, we translate pivot conditions for reactant and reaction matri- 
ces into the language of reaction networks. Note that reactant matrices M are non-negative 
matrices, so their pivots are necessarily positive pivots. 

Proposition 4.4. Let N be a square chemical reaction network with reactant matrix M and 
reaction matrix R. Then: 

1. Mij is a row pivot if and only if the i-th reactant contains only species j (in other words, 
the i-th reaction is kXj — ?• • ). 

2. Mij is a column pivot if and only if species j appears in the reactant complex of the i-th 
reaction and in no other reactant complex. 

3. Rij is a row pivot if and only if only species j is changed in reaction i. 

4. Rij is a column pivot if and only if species j is changed only in reaction i. 
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5. Both Mij and Rij are row pivots if and only if the i-th reaction is of the form 

ajXj — )• bjXj where ^ aj bj 

Moreover, Rij is a positive pivot if aj > bj and a negative pivot if aj < bj. If Mij and 
Rij are positive row pivots then we will say that i is a pivotal reaction. 

6. Both Mij and Rij are column pivots if and only if species Xi appears only in reaction j . 
Moreover, Rij is a negative column pivot if and only if Xi is a self-catalyst in reaction j. 

7. Suppose that for some i,j, both Mij and Rij are positive pivots. Then 

Or(iV) = sign(det(M) det(i?)) = sign{det{M^^) det{R^^)) . 
If Mij is a (positive) pivot and Rij is a negative pivot, then 

Or{N) = sign(det(M)det(i?)) = - sign(det(M*'^) det(i?*'^)) . 

Proof. This proposition follows easily from definitions and elementary linear algebra. □ 

4.2 Analysis of embedded networks 

Theorem 14.91 in this section will tell us that instead of analyzing all s-square subnetworks Gi of 
a network 0, we can remove the outflow reactions from Gi and then remove the outflow species 
from the remaining reactions, and analyze those simpler 'embedded networks.' First we must 
define what is meant by removing a reaction or a species. 

Definition 4.5. Let (3 = {S,C,TZ} be a chemical reaction network. 

1. Consider a subset of the species S C S, a subset of the complexes G C C, and a subset of 
the reactions R dlZ. 

• The restriction of R to S, denoted by R\s, is the multiset of reactions obtained by 
taking the reactions in R and removing all species not in S. 

• The restriction of C to R, denoted by G\r, is the set of (reactant and product) 
complexes of the reactions in R. 

• The restriction of S to C , denoted by S\c, is the set of species that are in the 
complexes in G . 

2. The network obtained from (5 by removing a subset of species {Xi} C S is the network 

{5\{XJ, C\^\,,,,^^, • 

3. The network obtained from (5 by removing a set of reactions {y —?■ y'} <ZlZ is the subnet- 
work 

{S\c\^\^^^^,y C\n\{y^y'}, n\{y^y'}^ . 
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Note that -RI5 is viewed as a multiset (a generalized set in which an element may appear more 
than once), so that in the following definition, embedded networks are indeed square. 

Definition 4.6. Let (5 = {S,C,TZ} be a chemical reaction network with s species, and let G 
be its non-flow subnetwork. For k = 0, 1, . . . ,s, a /c-square embedded network of & (or of 
G), which is defined by two sets, a k-subset of the species, S = {Xj^,Xj2, • • • i^i^} C S, and a 
k-subset of the non-flow reactions, R = {Rj^, Rj^-, ■ ■ ■ , Rju) ^ ^> involve all species of S, 
is the k-square network {S,C\ji^^, R\s) consisting of the k reactions R\s- 

The 1-square embedded networks with negative orientation are precisely the reactions of the 
form aXj bXj with 1 < a < 6. Therefore, a network contains a negatively oriented 1-square 
embedded reaction if and only if it contains a self-catalyzing reaction. 

The following proposition will be used in the proofs of Theorems 14.91 and 14.111 

Proposition 4.7. Let N be a square network that contains one of the following: 

1. an outflow reaction Xj — ?• or a generalized outflow reaction aXi bXi (where < 6 < 
a), or 

2. a species that appears in exactly one reaction, and the reaction is not self- catalyzing. 

Then Oi{N) = 0, or the network N'^^'^ obtained from N by removing the aforementioned reaction 
and species is a square network and Or{N) = Or{N^^'^). 

Proof. Part 1 follows easily from Proposition 14.41 &s such a reaction is a pivotal reaction. Now, 
part 2 is the case when some species Xj appears in exactly one reaction y y' in N. If Xj is 
in the product complex y' , then Or(iV) = by Lemma 13.101 In the remaining case {Xj is in 
the reactant complex y), then y ^ y' is a pivotal reaction, so Or(A^) = Or(A^'''''^). □ 

The following lemma generalizes what we saw in Example 14.11 in the CFSTR setting: each 
s-square subnetwork Gi of 6 defines an embedded network Ni of G, and they have the same 
orientation. 

Lemma 4.8. Let & be a chemical reaction network with s species, and let G denote the sub- 
network of non-flow reactions. Consider the map from s-square subnetworks Gi of <S that do 
not contain inflow reactions to square embedded networks of G, which maps Gi to the network 
Ni obtained from Gi by first removing all outflow reactions Xi^ — > (for I = 1,2, . . . , s — k) and 
then removing those outflow species Xi^ from the remaining non-flow reactions of Gi. Then this 
map is surjective, and 

Or(G,) = Or(iVi) . (4) 
(If Ni is empty, its orientation is defined to be -hi.) 

Proof. The surjection is clear: both objects are defined by a subset of S and a subset of TZg 
whose cardinalities sum to s. Equation follows from noting that the outflow reactions are 
pivotal reactions, and then applying Proposition 14.71 □ 
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From Lemma 14.81 we obtain the following theorem, which wih be used in Step 111 in Algo- 
rithm 16.11 



Theorem 4.9. For a CFSTR <3, with G its subnetwork of non-flow reactions, the following 
are equivalent: 

1. & passes the Jacobian Criterion. 

2. Each square embedded network of G has non-negative orientation. 

3. G has no self-catalyzing reaction, and each square embedded network NofG which satisfies 
the following properties has non-negative orientation: 

• N has no outflows, inflows, generalized inflow reactions ^-OiXi, or generalized 
outflows aXi —7- bXi (where < b < a). 

• Each species appears in at least two reactions (where a pair of reversible reactions is 
counted only once). 

Proof. Parts 1 and 2 are equivalent by Lemma [3.4l and Lemma [4.8i Clearly, part 2 implies part 3, 
so it remains only to show the converse. It suffices to show that if is a square embedded 
network which violates one of the conditions of part 3, then either it has zero orientation or we 
can find a 'smaller' embedded network which satisfies part 3 and has the same orientation as 
N. In the case of infiows and generalized inflows, the orientation of N is zero. All remaining 
cases are handled by Proposition 14.71 □ 

The following corollary follows easily from Theorem 14.91 

Corollary 4.10. If a CFSTR contains a self-catalyzing reaction, then it fails the Jacobian 
Criterion. Conversely, if a CFSTR does not contain a self-catalyzing reaction, then the only 
square embedded networks in Theorem \4-9\ that must be checked are those of size at least two. 

Corollary 14.101 allows us to see quickly by inspection that the network introduced in Example l2.4l 
(and analyzed in Examples 13.21 13.71 and 14. ip fails the Jacobian Criterion. 

4.3 A network reduction theorem 

This subsection concerns those reactions or species that can be removed from a CFSTR without 
affecting whether the network passes the Jacobian Criterion. The following theorem will be the 
foundation of Step 11 of Algorithm 16.11 

Theorem 4.11. Let & be a CFSTR, and let 0''°'^ be a network obtained from & by performing 
any of the following operations: 

1. Remove any inflow or outflow reactions. 

2. Remove any generalized inflow reactions — )• ajXj. 



13 



3. Remove reactions containing only one species aXi — t- hXi for which < b < a. 

4. Remove any species that is not a reactant species. 

5. Remove any species that appears in exactly one reaction (or a pair of reversible reactions), 
unless the species is a self-catalyst of the reaction. 

Then (5 passes the Jacobian Criterion if and only if all square embedded networks of (3^'^'^ have 
nonnegative orientation. 

Proof. Let G (respectively, C^'i) denote the subnetwork of non-flow reactions of & (respectively, 
0rcd^_ By Theorem 14.91 we need only show that all square embedded networks of & have non- 
negative orientation if and only if all square embedded networks of G^^'^ do. First, we consider 
the operations 1 through 3, in which is obtained from 05 by removing a reaction y y' . 
In this case, 

{SENS of G} = { SENS of G''"^} U { SENs of G that involve y ^ y'} , (5) 

where 'SEN' is shorthand for 'square embedded network.' So, we need only show that if Ni is 
a square embedded network of G, then either Ni has zero orientation, or there exists a square 
embedded network N^^^ of G^'^'^ which has the same orientation as Ni. In the case that the 
reaction y ^ y' is an inflow or a generalized inflow, Or(A^j) = 0. So, the remaining cases are 
when y — 7- is an outflow reaction Xj ^ or is of the form aXj — )■ bXj (where 1 < b < a); 
this is handled by Proposition 14.71 

Now, only parts 4 and 5 remain to be proven. Here, ^5''^'^ is obtained from (3 by removing 
a species Xj. The analogue to ([5]) now is that the set of square embedded networks of C'^'^ are 
precisely the square embedded networks of G that do not involve species Xj. Again, we must 
show that a square embedded network Ni of G either has zero orientation, or there exists a 
square embedded network N^^^ of G^^'^ with the same orientation as Ni. If Xj is not a reactant 
species of G (in the case of part 4), then Ni has zero orientation by Lemma [3. 101 Finally, if Xj 
satisfies the condition in part 5, then Xj appears in exactly one reaction y y' in Ni. This 
case is completed by applying Proposition 14.71 This completes the proof. □ 

5 Networks with total molecularity at most two do not admit 
multiple steady states 

In this section, we prove the following theorem, which concerns networks in which the maximum 
total molecularity (Definition 12. 2p is at most two. 

Theorem 5.1. If all species of a CFSTR have total molecularity less than or equal to two, then 
it passes the Jacobian Criterion, and therefore does not admit multiple steady states. 

Theorem 15.11 has the following corollary. 
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Corollary 5.2. For any weakly -reversible chemical reaction network (not necessarily a CFSTR) 
in which all species have total molecularity less than or equal to two, the network does not admit 
multiple steady states. 

Remark 5.3. Theorem \5.1\ cannot be improved to higher total molecularity. For example, a 
CFSTR that contains the reaction A — > 2A fails the Jacobian Criterion (by Corollary \4. 10^ . 

Remark 5.4. Theorem \5.1\ can be proven by applying the Species- Reaction Craph Theorem of 
Craciun and Feinberg Q Theorem 1.1]. Rather than introducing the terminology and results of 
species-reaction graphs, we give a direct proof here. 

Our proof of Theorem 15.11 will proceed by showing that for a chemical reaction network 
that has total molecularity at most two, if a square embedded network has nonzero orientation, 
then it in fact has positive orientation. The proof requires the following lemmata, which will 
be used to reduce to studying square embedded networks, and then to elucidate what values 
can be taken by the determinants of the reactant and reaction matrices of such networks. 

Lemma 5.5. Let & be a chemical reaction network whose species all have total molecularity 
less than or equal to some m € Z>o. Then the same holds for any embedded network of &. 

The proof of Lemma 15.51 follows immediately from the definitions of total molecularity and 
embedded networks. 

Lemma 5.6. Let N be a coupled s-square network with at least two reactions (none of which 
are inflows or outflows) such that all species have total molecularity at most two. Assume that 
N has nonzero orientation. Then each reaction of N has the form of A-\- B or A B. 

Proof. We first claim that if is a coupled, s-square network in which all species have total 
molecularity two, then each reaction of must contain exactly two molecules. This is because 
no reaction can have only one molecule (because inflows and outflows are not allowed). So, each 
of the s reactions must have exactly two molecules, and it follows that the total molecularity of 
each species must be two. Now, a reaction of the form ^ A -\- B would yield zero orientation 
for the network (because the corresponding row of the reactant matrix would be zero), and a 
reaction 2A — )■ is not allowed, because then the network would decouple into this reaction 
and the remaining subnetwork that does not contain A. So, only two types of reactions are 
possible: A + B ^ and A^ B. □ 

Lemma 5.7. Consider an n x n-matrix Q G { — 1, 0, 1}"-^"- such that 

1. all diagonal entries are +1, 

2. the entries above the diagonal (mod n), Qjj+i(mod n)> o,re or ±1, and 

3. all other entries are 0. 

Then dei{Q) > 0. (In fact, the determinant is 0, 1, or 2.) 



15 



Proof. The determinant is det(Q) = 1 ± ll'i=iiQi,i+iimod n)) ^ {0, 1, 2}. □ 

We now can prove Theorem 15.11 and Corollary 15. 2[ 

Proof of Theorem \5.1[ By Theorem 14. 9| Lemma 14.81 and Lemma 15.51 we need only show that 
if is a square network with nonzero orientation and with maximum total molecularity at 
most two, then det{Mi^) det(i?jv) ^ 0. Let s denote the number of species of N. Now, if N 
decouples, then we could reduce to analyzing the two subnetworks separately both of which 
have fewer than s species; the network N passes the Jacobian Criterion if and only if the two 
subnetworks do. The s = 1 case is trivial (only 2A — >■ is allowed, and in this case both reaction 
and reactant matrices are the same, so their determinants are the same), so we may proceed 
by induction on s in the case of a coupled network A^. 

Now Lemma 15.61 applies, so we know that the reactions of N are of the form A + B ^ 
and A ^ B. By relabelling the species, we can ensure that the first reaction contains species 
Xi and X2, the second reaction contains X2 and X3, and so on. Therefore, both the reactant 
matrix and the reaction matrix now have the form of the matrix in Lemma 15.71 By the 
lemma, both determinants are non-negative, so the orientation of the network is either zero 
or positive. 

Finally, Lemma l3 . 41 implies that a CFSTR network that passes the Jacobian Criterion does 
not admit multiple steady states. □ 

Proof of Corollary \5.B . For a weakly-reversible network G that has total molecularity at most 
two, we consider the CFSTR G obtained by augmenting the network G by all outflows. The- 
orem 15.11 implies that G passes the Jacobian Criterion. Finally, results due to Craciun and 
Feinberg allow us to conclude that the such a network precludes multiple steady states [HI 
Theorem 8.2 and Corollary 8.3]. □ 



6 Procedure for simplifying the Jacobian Criterion 

Algorithm 16. 11 is the main result of this article. It combines the results in previous sections into 
a procedure that simplifies the implementation of the Jacobian Criterion. After pre-processing 
in Step I, the algorithm reduces a chemical reaction network to a simpler network in Step II, 
and then decreases the number of embedded networks of the reduced network that must be 
examined in Step III. 

Algorithm 6.1 (Simplifying the Jacobian Criterion). 
Step I (Pre-processing step). 
Input: a CFSTR &. 

Output: either whether & passes the Jacobian Criterion, or & (to be passed to Step II). 

1. If 6 has a self- catalyzing reaction (an embedded reaction aXj — t- bXj with 1 < a < b), 
then (3 automatically fails the Jacobian Criterion. 

2. If each species of © has total molecularity at most two, then (S automatically passes the 
Jacobian Criterion. 
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If neither applies, then pass <3 to Step II. 

Step II (Reduction step). 

Input: a CFSTR © that passes Step I. 

Output: a non-flow chemical reaction network Gred such that the CFSTR G^ed passes the 
Jacobian Criterion if and only if <3 does. 

Perform the following operations on the network (S, in any order, as long as any one can be 
completed. 



1. 


Remove 


any inflow or outflow reaction. 


2. 


Remove 


any generalized inflow reaction — )• ^jOjXj. 


3. 


Remove 


any reaction that contains only one species. 


4- 


Remove 


any species that is not a reactant species. 


5. 


Remove any species that appears in only one reaction (or in a pair of reversible reactions 
and in no other reactions). In particular, remove any species whose total molecularity is 
one. 


6. 


Remove 


one copy of any repeated reactions. 



If the reduced network is the empty network or all species in the reduced network have total 
molecularity at most two, then passes the Jacobian Criterion. Otherwise, pass the reduced 
network G^'^^ to Step III. 

Step III (Analysis of embedded networks step). 

Input: a chemical reaction network G (such as the output of Step II) which passes Step I. 
Output: whether G passes the Jacobian Criterion (if G is the output of Step II, then this step 
will determine whether the input © of Step II passes the Jacobian Criterion). 
Consider all coupled square embedded networks N of G with the following properties: 

1. N contains none of the following reactions: 

• inflow or outflow reactions, 

• generalized inflow reactions — ?• ^^aiXj, 

• reactions containing only one species, 

• both a reaction and its reverse reaction, 

• two reactions that share the same reactant complex. 

2. N has at least two species. 

3. Each species of N is a reactant species. 

4. Each species appears in at least two reactions (where pairs of reversible reactions are 
counted only once). 
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5. At least one species has total molecularity three or more. 

6. N has nonzero orientation. 

Output: The network G passes the Jacobian Criterion if and only if all such square embedded 
networks N have nonnegative orientation. 

Proof of correctness of Alaorithm \6.1\ Part 1 of Step I follows from Corollary 14.101 part 
2 follows from Theorem I5.1i Step II is valid by Theorem 14.111 Parts 1 through 4 of Step III 
follow from Lemma 13.1^ Corollary 14.101 3.nd Theorem 14.91 Part 5 of the Step III follows from 
Theorem 15.11 and part 6 simply avoids embedded networks with zero orientation. □ 

Remark 6.2. Instead of applying Step III to a reduced network that is the output of Step II, 
one can apply other tests to check the Jacobian Criterion. For instance, one can apply the 
Species- Reaction Graph Criterion to the reduced network or directly compute the expansion 
of the Jacobian determinant using available software /731 [HJ The reduced network typically 
contains fewer species and fewer reactions than the original network, so computer computation 
will be more tractable. 

Remark 6.3. In practice, the number of square embedded networks that must be examined in 
Step III of Algorithm 16. il is relatively small compared to the number of s-square subnetworks 
Gi in Lemma \3.4\ Comparisons of the number of square subnetworks Gi versus the number 
of square embedded networks to be examined appear in the last two columns of Table [7] and in 
the comments of Table In addition, the determinants of such square embedded networks ( as 
given in Lemma \3.3\) that must be computed are of smaller matrices than those of Gi. 

Remark 6.4. We propose that Alaorithm \6.1\ could aid in enumerating chemical reaction net- 
works that admit multiple steady states. Recall that failing the Jacobian Criterion is a necessary 
condition for admitting multiple steady states (Lemma \3.4% so let us consider the problem of 
identifying the set of networks that fail the Jacobian Criterion. More specifically, we now define 
a bimolecular network to be a network in which all complexes contain at most two molecules, 
and then pose the following question: which bimolecular CFSTRs fail the Jacobian Criterion? 

Accordingly, we enumerated all bimolecular CFSTRs that contain two pairs of reversible 
non-flow reactions and we found that there are 386 such networks. Of these 386 networks, 14-2 
have a maximum total molecularity of two or fewer, and thus immediately pass the Jacobian 
Criterion. It is interesting to note that for these 1^2 networks, several distinct species typically 
have the same total molecularity. Networks that have such symmetry are computationally the 
most expensive to enumerate Jllf . As we are interested in listing only those networks that fail 
the Jacobian Criterion, we can avoid these computationally expensive networks. 

Moreover, of the remaining 244 networks, applying Step II of the Alaorithm \6.1\ immediately 
reduces the number of networks that need to be examined. In other words, many networks are 
related in whether or not they pass the Jacobian Criterion and thus studying all 244 networks 
is redundant. When we discard networks in which some species has total molecularity equal to 
one, we are left with only 55 networks whose multistationarity needs to be investigated. Details 
of our enumeration of bimolecular CFSTRs appear in 120^ . Note that a similar enumeration 
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of small himolecular reaction networks was undertaken by Deckard, Bergmann, and Sauro \1 jj/ . 
from which Pantea and Craciun sampled networks to compute the fraction of such networks that 
pass the Jacobian Criterion \2'A Figure 1]. 



7 Examples 

The examples in this section iUustrate that Algorithm [6TT] in some cases can determine whether a 
network passes the Jacobian Criterion more quickly than applying the Species- Reaction Graph 
Criterion of Craciun and Feinberg [7j, Deficiency Theory [14^ I15j. or the results of Banaji and 
Craciun [2]. Moreover, we are able to analyze examples in [51[71 lim i25|. which appear in Tables [T] 
and[2j In addition, we demonstrate that instead of checking the orientations of all the s-square 
networks Gi from Lemma 13.41 (which entails computing determinants of many s x s-matrices), 
we need only compute the orientations of certain square embedded networks, which are both 
smaller and fewer in number. The final columns in both Tables [T] and [2] illustrate the great 
reduction in the number of subnetworks that must be examined. We begin with examples in 
which no orientations must be computed, as checking the Jacobian Criterion sometimes can be 
performed by inspection. 



Example 7.1 (Networks containing a self-catalyzing reaction). Recall that Corollary \4-10 states 



that if a chemical reaction network contains a self-catalyzing reaction, then it fails the Jacobian 



Criterion. One such network is the one examined in Examples 3.1 and \J^.1\ Other networks 
that are resolved easily by Step I of Algorithm \6.1\ include networks (vi), (vii), and (viii) in 
Table □ 

In addition to the self-catalyzing networks, the other networks that are resolved in Step I 
of Algorithm 16.11 are those with total molecularity at most two. 

Example 7.2 (Networks with total molecularity at most two). Consider the following network 
examined by Gnacadja et al. 

L + R^LRA + R^AR (6) 
L + T ^ LT A + T ^ AT . 

This ligand-receptor-antagonist-trap network summarizes the four pairs of binding and dissociat- 
ing reactions that occur in the presence of two drugs, a 'receptor antagonist' and a 'trap, ' which 
are used to treat patients with rheumatoid arthritis. In their work, Gnacadja et al. apply Defi- 
ciency Theory to determine the existence and uniqueness of a steady state IT1\ Proposition 2]. 
If one is interested only in precluding multiple steady states, then the approach advocated here 
can be performed more quickly. It is easily seen that the total molecularity of each of the eight 
species, L, R, A, T, LR, AR, LT, and AT, is no more than two. Therefore, Theorem I5.il 
implies that a CFSTR version of this network passes the Jacobian Criterion, and does not ad- 
mit multiple steady states. Moreover, as the underlying network ([6]) is weakly-reversible, we can 
conclude that multiple steady states is also ruled out in the non-GFSTR setting J^. 
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C^A 






see Ex. [TJ 




(v) 


A + B^F 














A + C^G 


No 


No 


See Ex. 1731 


(8+7) = 6435 
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C + D^B 












C + E^D 












(vi) 


A + B^2A 


No 


Yes 


Self-catalyzing, 
see Ex. [Q 


m = 6 




(vii) 


2A + Bt^3A 


No 


No 


Self-catalyzing, 
see Ex. O 


(^f ) = 6 




(viii) 


A + 2B '=^3A 


No 


Yes 


Self-catalyzing, 
see Ex. O 


CV) = 6 





Table 1: Networks examined in Table 1 of Schlosser and Feinberg [25]. Listed are the non-flow 
networks G, whether the CFSTR network G passes the Jacobian Criterion (JC), whether G 
precludes multiple steady states (MSS), additional comments, the number of s-square networks 
Gi of G that do not contain inflow reactions (here, Tq = s; see Lemma l3.4p . and the number 
of square embedded networks (SENs) which must be considered in Algorithm 16.11 The entry 
in the last column is left blank (— ) for networks that can be determined to pass or to fail 
the Jacobian Criterion due to a self-catalyzing reaction or low total molecularity. Note that 
the number of such SENs is much fewer than the number of subnetworks Gi, resulting in a 
substantial reduction in the number and size of the orientations that must be computed. Many 
of these networks were also analyzed by Craciun and Feinberg [5l Table 1.1], Banaji, Donnell, 
and Baigent [3] and Craciun, Helton, and Williams [9]. 
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G 


G does 




Network 


(G) 


passes 
JC 


not admit 
MSS 


Remarks 


(i) 


E + S^ES^E + P 


Yes 


Yes 


Total molecular ity 
< 2, see Ex. O 


(ii) 


E + S^ES^E + P 
E + I^EI 


Yes 


Yes 


Reduced network has 
TM < 2, see Ex. O 


(iii) 


E + S^ES^E + P 
ES + I^ ESI 


Yes 


Yes 


Reduced network has 
TM < 2, see Ex. O 


(iv) 


E + S^ES^E + P 






Only 9 SENs 




E + I^EI 


No 


No 


to examine, 




ES + I^ ESI ^EI + S 






see Example 17.51 


(v) 


E + SiU ESi 

S2 + ESi ^ ES1S2 ^E + P 


Yes 


Yes 


Total molecularity 
< 2, see Ex. [7J 


(vi) 


E + Si^ESi E + S2^ ES2 






Only 14 SENs 




S2 + ESi ^ ES1S2 ^ iSi + ES2 


No 


No 


to examine, 




ES1S2 ^E + P 






see Example 17.51 


(vii) 


Si + El ^ EiSi 






Only 1 SEN 




S2 ~t~ El Si E1S1S2 — ^ Pi ~t~ El 


Yes 


Yes 


to examine. 




E2 + 5*2 ^ E2S2 — ^ 25"! + E2 






see Example 17.51 


(viii) 


5*1 + El ^ El Si 

S2 + EiSi ^ E1S1S2 — >■ i^i + i?i 






Only 1 SEN 




5*2 + -E2 ^ E2S2 


No 


No 


to examine. 




5*3 + E2S2 ^ E2S2S3 — 7- P2 + -E'2 






see Example 17.51 




53 + E3 ^ E353 ^ 25i + ^3 








(ix) 


5*1 + i?i ^ El 5*1 

5*2 + Ei^i ^ E1S1S2 — )• Pi + El 

5*2 + E^2 E2S2 






Only 1 SEN 




S3 + E2S2 ^ E2S2S3 — 7- P2 + E2 


Yes 


Yes 


to examine. 




5*3 + £^3 ^ E;3'S'3 






see Example 17.51 




54 + £^35'3 ^ £^3S'3S'4 ^ P3 + E3 










5*4 + E/4 ^ E4S4 — 25"! + E'4 









Table 2: The enzyme catalysis networks examined in Table 1 of Craciun, Tang, and Fein- 
berg [To]. Listed are the underlying networks G, whether the CFSTR network G passes the 
Jacobian Criterion (JC), whether G precludes multiple steady states (MSS), and additional 
comments, including the number of square embedded networks (SENs) that need to be exam- 
ined for those networks for which total molecularity (TM) is not sufficient for verifying the 
Jacobian Criterion. 
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Other networks with total molecularity at most two include the network (iv) in Table{l\ and 
the enzyme catalysis networks (i) and (iv) in Tahle\^ All of these networks can he seen by 
inspection to pass the Jacobian Criterion, without appealing to the SR Graph Criterion. 

The networks in Example 17.21 had total molecularity at most two, and so could be seen 
immediately to pass the Jacobian Criterion. The networks in the next example do not have 
this property until after applying Step II of Algorithm 16. 1[ 

Example 7.3 (Networks with total molecularity at most two after reduction). Consider the 
CFSTR introduced by Banaji and Craciun whose non-flow subnetwork is: 

A + B + C A + B + C A + B (7) 

Banaji and Craciun showed that the CFSTR fails the Species-Reaction Graph Criterion, but 
nevertheless passes the Jacobian Criterion, because the stoichiometric matrix is 'strongly sign 
determined' |^ ^6-2]. We now show that Step II of Alaorithm \6.1] can easily be performed on 
this network, to conclude that it passes the Jacobian Criterion. 

In Step I of Algorithm \6.1[ we note that the total molecularity of species A and B in net- 
work ([7]) are both three, and no self-catalyzing reactions exist. So, we begin Step II by removing 
species D, E, and F, each of which has total molecularity equal to one. The resulting network 
is: 

Q^A + B + C Q^A + B + C Qt^A + B 

Now the first two pairs of reversible reactions are the same, so we remove one copy. Additionally, 
we remove the remaining two generalized inflow reactions, to obtain: 

A + B + C A + B 

Clearly, each species in this reduced network has total molecularity no more than two, so we 
conclude that the original CFSTR passes the Jacobian Criterion. Two other networks whose 
reduced networks can be obtained easily and have total molecularity at most two are the enzyme 
catalysis networks with inhibition (ii) and (Hi) in Table\^ 

While the above examples did not need to proceed beyond Step II of Algorithm 16.11 the 
networks in the following example do require this step. We will see that the algorithm greatly 
reduces the number of subnetworks that must be checked. 

Example 7.4 (Networks requiring the full algorithm). Consider the reversible network (i) in 
Table\^ its non-flow subnetwork is: 

A + B^P B + C^Q C^2A (8) 

The output of Step II of Alaorithm \6.1\ is the following network: 

A + B^Q B + C ^0 C^2A 
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As only species A has total molecularity greater than two, each square embedded network we 
consider must contain both A + B and one of the reversible reactions C ^ 2A (because 
we need not consider networks with both directions of a reversible reaction). However, such a 
2-square embedded network would have a species (B or C) that appears in only one complex. 
Therefore, we need only consider the two 3-square embedded networks, which are the following: 

{A + B^O B + C ^0 C ^ 2A} , and (9) 
{A + B^O B + C^O C^2A} (10) 

The first embedded network Q has negative orientation, and the second network (jlOp has pos- 
itive orientation, so a CFSTR version of the original network ([8|) fails the Jacobian Criterion. 
Note that our algorithm reduced the Jacobian Criterion from examining the (^+^) = 462 s- 
square subnetworks (the Gi in Lemma \3.4^ to two embedded networks. Networks (ii), (Hi), 
and (v) in Table [1\ can be analyzed similarly: they require examining only two, two, and three 
square embedded networks, respectively. 

Our final example network requires more complicated analysis, but the number of square 
embedded networks that must be checked is nevertheless small. 

Example 7.5. Consider network (iv) in Table\M After removing species 'P' in Step II and 
then relabeling the six remaining species, the subnetwork of non-flow reactions is: 

A + B^C C^A A + E^F C + E^G G^F + B 

Only one reaction, G ^ A, is not reversible; we label the remaining four reversible reactions 
as Ri through R4. We now show that instead of checking that the (^y) = 11440 7 -square 
subnetworks of the original CFSTR (with 7 species and 9 non-flow reactions) have non-negative 
orientation, Algorithm \6.1\ allows us to examine only nine square embedded networks. 

Let N denote a square embedded network that must be considered according to Step III 
of Algorithm The total molecularities of the species {A,B,G,E,F,G) are respectively 

(3,2,3,2,2,2). So species A together with all its reactions (Ri, R2, R3) or species C together 
with all its reactions (Ri, R2, Ra) must be in N. In particular, both reactions Ri and R2 must 
appear in N, and both A and G must take part so that R2 is not a flow reaction. Next, in order 
for species C not to appear twice as a reactant, Ri must appear in N in the forward direction 
A + B^G. 

Now, the rows corresponding to reactions Ri and R2 in the reaction matrix of N must be 
linearly independent, so the species B must be in N. Thus, we also must include reaction R^ 
so that B has total molecularity greater than 1. Finally, so that R^ is not a flow reaction, we 
must include species F or G. 

Therefore, N must be a 4-square or a 5-square embedded network, and the only possible 
4-square networks are the following: 

{A + B^G G^A B^G G^Cjand 
{A + B^G G^A F + B^O A^ F} , 
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where the second network represents two square embedded networks. Additionally, the only 
possible 5-square networks are 

{A + B^C C^A A^F G^C F + B ^ G} , 

{A + B C ^ A A + E ^0 G + E^G G ^ B} , and 

{A + B^G G^A A + E^F G + E^O F + B ^ 0}. 

Here, the second network represents three networks (the network in which G appears twice as 
a reactant is not permitted), and the third network represents two networks. Therefore, there 
are 9 possible square embedded networks N . Using similar analysis, network (v) of Table{l\ and 
networks (vi), (vii), (viii), and (ix) of TablelM can be shown to require checking only 3, 14, 1, 
1, and 1 square embedded networks, respectively. 
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